Multiscale quantum criticality: Pomeranchuk instability in isotropic metals 
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As a paradigmatic example of multiscale quantum criticality, we consider the Pomeranchuk in- 
stability of an isotropic Fermi liquid in two spatial dimensions, d = 2. The corresponding Ginzburg- 
Landau theory for the quadrupolar fluctuations of the Fermi surface consists of two coupled modes, 
critical at the same point, and characterized by different dynamical exponents: one being ballistic 
with dynamical exponent z = 2 and the other one is Landau damped with z = 3, thus giving rise 
to multiple dynamical scales. We find that at temperature T — 0, the ballistic mode governs the 
low-energy structure of the theory as it possesses the smaller effective dimension d + z. Its self- 
interaction leads to logarithmic singularities, which we treat with the help of the renormalization 
group. At finite temperature, the coexistence of two different dynamical scales gives rise to a mod- 
ified quantum-to-classical crossover. It extends over a parametrically large regime with intricate 
interactions of quantum and classical fluctuations leading to a universal T dependence of the corre- 
lation length independent of the interaction amplitude. The multiple scales are also reflected in the 
phase diagram and in the critical thermodynamics. In particular, we find that the latter cannot be 
interpreted in terms of only a single dynamical exponent whereas, e.g., the critical specific heat is 
determined by the z = 3 mode, the critical compressibility is found to be dominated by the z = 2 
fluctuations. 



PACS numbers: 



I. INTRODUCTION 



At a quantum phase transition, statics and dynamics 
of critical fluctuations are inseparably entangled. This 
is manifest in the dependence of static and thermody- 
namic properties on the dynamical exponent z, that al- 
lows to associate a thermal length, £r ~ T~ x / Z , with 
temperature T. This thermal length is at the origin of 
the quantum-to-classical crossover where critical fluctua- 
tions change their character^ Whereas critical bosonic 
modes with momenta higher than the inverse thermal 
length are of a quantum-mechanical nature, modes with 
smaller momenta behave effectively classical. This cross- 
over also reveals itself in the phase diagram spanned by 
the control parameter of the quantum phase transition, 
r, and temperature, T. The behavior of, e.g., thermody- 
namic properties changes qualitatively when the correla- 
tion length £ is of the same order as the thermal length 

The identification of the critical fluctuations and their 
dynamics close to a putative quantum phase transition 
in a material is a nontrivial task. Frequently low-energy 
fluctuations characterized by different dynamics coexist 
and interact. In fact, this is generally the case in crit- 
ical metals where the dynamics of the order parameter 
coexists with the dynamics of the fermionic quasiparti- 
cles. The standard Hertz-Millis-Moriya model^££ for 
a magnetic instability in a metal assumes that only the 
dynamics of the paramagnons is important for the char- 
acterization of critical properties. Although this assump- 
tion has proved to be successful in interpreting a number 
of experiments, it has failed to explain consistently the 
properties of the heavy fermion compounds CeCu6- x Au x 



and YbRh2Si2^^ In particular, it has been argued that 
the latter system is characterized by the presence of mul- 
tiple scales^ that might be due to the presence of coex- 
isting dynamics. Theoretically, it is now established that 
the assumptions of the Hertz-Millis-Moriya approach fail 
in the case of the ferromagnetic instability— The coupling 
of the ferromagnetic fluctuations to the fermionic degrees 
of freedom gives rise to singular corrections that prevent 
the formulation of an effective theory exclusively in terms 
of order-parameter fluctuations i^ ' 10 ' 11 ' 12 The presence of 
multiple dynamics is also important for the antifcrromag- 
netic quantum phase transition in metals in spatial di- 
mension d = 2; the fermions mediate an effectively long- 
ranged interaction among the paramagnons resulting in a 
strongly coupled theory with multiple dynamical scales^ 

Although quantum criticality with multiple dynamical 
scales seems to be rather common, systematic investiga- 
tions are rare and the generic features of this problem 
are therefore only poorly understood. The presence of 
coexisting dynamics and, consequently, different dynam- 
ical exponents complicates the analysis of experimental 
data considerably, as simple scaling relations are bound 
to fail. Generally speaking, different physical quantities 
might be dominated by one or the other dynamics such 
that the interpretation in terms of a single exponent z 
results in apparent inconsistencies. On the theoretical 
side, a critical field theory with multiple exponents z 
promises to exhibit a rich structure. The modes with 
the larger dynamical exponent z> have a larger phase 
space available and are thus expected to dominate, e.g., 
the specific heat. The modes with the smaller exponent 
z<, on the other hand, have a smaller effective dimension 
d+z < and might therefore trigger infrared singularities in 
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perturbative loop corrections. In addition, multiple dy- 
namical exponents imply multiple thermal scales £r with 
the concomitant crossovers. Instead of a single quantum- 
to-classical crossover one expects an extended crossover 
regime with coexisting and interacting quantum and clas- 
sical fluctuations. 

In the present work, we theoretically investigate these 
phenomena in a simple model exhibiting multiscale quan- 
tum criticality. Our aim here is not to go beyond the 
Hertz-Millis-Moriya framework. Instead, we still limit 
ourselves to an effective bosonic theory where the fluctu- 
ations are however characterized by multiple dynamical 
scales. To be specific, we consider the effective Ginzburg- 
Landau theory for the Pomcranchuk instability^ in an 
isotropic metal as proposed by Oganesyan et a/.— The or- 
der parameter is given by the shear modes of the Fermi 
sphere that drive a spontaneous deformation of the Fermi 
surface at the instability resulting in a "nematic" state. 
In d = 2, there are two critical bosonic shear modes char- 
acterized by different dynamical exponents. Whereas one 
is damped by particle-hole excitations in the metal and 
has a dynamical exponent z > = 3, the other is ballistic 
and undamped with z< = 2. We confine ourselves to two 
spatial dimension, d = 2, where symmetry consideration 
allow the Pomeranchuk transition to be of second order 
and fluctuation corrections are pronounced. 

The Pomcranchuk instability and its criticality has re- 
ceived renewed interest in recent years. In analogy to 
paramagnons close to a magnetic quantum phase tran- 
sition in metals, the critical collective shear fluctua- 
tions close to a Pomeranchuk quantum critical point give 
rise to a singular interaction among the fermionic de- 
grees of freedom leading to a breakdown of Fermi-liquid 
theory i 15 i 16 i 17 The structure of the critical theory, how- 
ever, depends on the presence of an underlying crystal 
lattice. Whereas in an isotropic metal the Fermi-liquid 
properties are destroyed over the full Fermi surface^ 
well-defined quasiparticles partly survive if a crystal lat- 
tice breaks rotational symmetry from the outset i 16 i 17 
Moreover, the aforementioned undamped z < = 2 mode 
only becomes critical in the isotropic case but is always 
gapped if a crystal is present. A pecularity of the ne- 
matic, i.e, Pomeranchuk-ordcrcd phase in an isotropic 
medium is the existence of a Landau-damped Goldstonc 
mode that also leads to singular Fermi-liquid corrections 
thus providing a rare example of an extended non-Fermi- 
liquid phased An extensive investigation of the mean- 
field properties of the Pomcranchuk transition has been 
carried out by Yamase and collaborators i 18 i 19 ' 20 i 21 A non- 
perturbative analysis of Pomeranchuk quantum critical- 
ity was presented in Refs. I^I^EmI . and a discussion 
of the Pomeranchuk instability from the perspective of 
Fermi- liquid theory was given by one of the authors 3i 
Quintanilla et aL 26 i 27 presented a careful study of as to 
how the instability might arise from central interactions 
among quasiparticles. Finally, the influence of disorder 
on the transition was studied in Ref. [28l . Nematic or- 
dering was also discussed extensively in the context of 
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FIG. 1: Effective dimensions of the longitudinal and transver- 
sal shear fluctuations as a function of momentum q. The re- 
spective thermal momenta, ~ T 1 / 2 , separate two regimes 
where the fluctuations have either a quantum character, 
q > ^y 1 , with effective dimension d + z or a classical char- 
acter, q < ^y 1 , with effective dimension d — 2. There is an 
extended quantum-to-classical crossover (gray shaded), where 
the quantum regime of the transversal z< = 2 mode overlaps 
with the classical regime of the longitudinal mode z > = 3. 

cuprate superconductors, for a recent review see Ref. |29| . 
On the experimental side, a metallic nematic phase in 
Si'3Ru207 close to its metamagnetic transition has been 
reported by Borzi et ab 30 ^ 31 A Pomeranchuk instability 
in an isotropic metal, that we consider here, has not been 
realized so far but might be achievable in cold atom sys- 
tems close to a Feshbach resonance in a channel with 
higher angular momentum. 



A. Summary of main results 

In the literature, attention has been mainly focused 
on the damped z > — 3 component of the shear fluctua- 
tions. Oganesyan et alJ^ even argued that the undamped 
z < = 2 fluctuation mode plays no role in the critical the- 
ory because its dynamics seems to be irrelevant with re- 
spect to z> = 3 scaling; hence, the critical theory should 
be fully captured by a Gaussian theory as the effective 
dimension d + z> = 5 is larger than the upper critical 
dimension of the effective <f> 4 theory, d+ = 4. In contrast 
to this, we find that the scaling arguments of Ref. [T^ are 
misleading and that the z< = 2 mode is in fact instru- 
mental for the Pomeranchuk quantum phase transition 
in d — 2. The ballistic mode has an effective dimension 
d + z< = 4 = d£ , and it therefore generates important 
logarithmic singularities at T = in loop corrections. We 
treat these singularities with the help of the renormaliza- 
tion group (RG). The derived RG equations identify the 
universality class which bears signatures of the tensorial 
nature of the shear fluctuations and differs, in particular, 
from the ones with Ising as well as XY symmetry. By 
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solving the RG equation, we find, for example, that the 
correlation length £ diverges with vanishing distance to 
the quantum critical point, r — > + , according to 



T ' 1 



|T=0 



(log^) 4 / 9 ' 



(1) 



where A is some momentum cutoff. The characteristic 
logarithmic enhancement with exponent 4/9 explicitly 
demonstrates that interactions among the z < = 2 modes 
lead to qualitative corrections that are beyond Gaussian 
critical behavior. 

At any finite temperatures, the coexisting dynamics 
implies the presence of two thermal lengths, £^ ~ j 1 ~ 1 / z > 
and ^ ~ T _1 / 2< , that give rise to an interesting modi- 
fication of the quantum-to-classical crossover. Generally 
at criticality, r = 0, critical modes change their charac- 
ter at the momentum scale of the inverse thermal length, 
q ~ £f ■ For l ar g er momenta q > S,^ 1 , the fluctua- 
tions are essentially of a quantum-mechanical nature be- 
cause the associated Matsubara frequencies can be ap- 
proximated to be dense. For smaller momenta q < ^y 1 , 
on the other hand, fluctuations are effectively classical 
in the sense that they are dominated by the Matsubara 
zero mode. In the presence of only a single thermal length 
£r, an effective theory at criticality, r = 0, can be for- 
mulated that changes its character at the momentum 
scale q ~ from a quantum theory with dimension 
d + z to an effective classical theory with dimension d, 
which is also known as dimensional reductioni^^ In 
the present problem of the Pomeranchuk instability we 
find, however, that it is not possible to derive the crit- 
ical properties from an effective classical model so that 
the concept of dimensional reduction breaks down. As 
illustrated in Fig. [TJ the two thermal lengths give rise to 
an extended intermediate momentum regime where the 
z < = 2 fluctuation mode is still quantum and the z > = 3 
mode is already classical. Remarkably, it turns out that 
interactions between quantum and classical fluctuations 
with momenta belonging to this extended quantum-to- 
classical crossover regime determine the critical proper- 
ties at r — resulting in the inapplicability of dimen- 
sional reduction. At d + z < = 4 = d+, the z < = 2 
quantum fluctuations generate scale-dependent vertices 
in loop corrections, that are in turn probed by the clas- 
sical z> = 3 mode. This peculiar interplay of quantum 
and classical fluctuations yields, for example, a universal 
temperature dependence of the correlation length in the 
low-T limit, 



r 



lr=0 



CT, 



4/9' 



(2) 



where £ and T are measured in dimensionless units to 
be specified below. The T dependence is universal in 
the sense that it does not depend on the value of the 
$ 4 -interaction amplitude although it is triggered by the 
very same interaction. Note that the correlation length 
[Eq. (J3J)] is of the same order as the thermal length 
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FIG. 2: Phase diagram for the Pomeranchuk quantum phase 
transition in d = 2. The correlation length £ and the two ther- 
mal lengths, ^ ~ T~ 1/3 and ^ ~ T~ 1/2 , distinguish three 
regimes: (I) a Fermi-liquid regime where £ is the smallest 
scale, (II) an overlap regime where £ is sandwiched between 
the two thermal lengths, and (III) a quantum critical regime 
where £ ~ £y is the largest scale. The line labeled Tg{v) 
represents the Ginzburg temperature. 



T- 1 / 2 which is a clear manifestation for the break- 



down of dimensional reduction. 

The coexistence of two dynamical scales is also re- 
flected in the phase diagram and in the critical ther- 
modynamics. The two thermal lengths £^ < £J? divide 
the disordered region of the phase diagrams into three 
regimes, see Fig. [2] (I) a low-temperature regime for 
£ < £^ , (II) an overlap regime, £^ < £ < £yJ, and 
(III) a quantum critical regime, £^ ~ £. Neglecting loga- 
rithmic corrections, the two crossovers are located in the 
phase diagram at T ~ r 3 / 2 and T ~ r, respectively. In 
the present work, we limit ourselves to the disordered 
side of the phase diagram on the right-hand side of the 
Ginzburg temperature, Tq(t) (non-shaded in Fig. 

It is a priori not clear which of the two crossover scales 
is more important for thermodynamics. Our calculations 
reveal that it depends in fact on the particular thermo- 
dynamic quantity. Due to their larger phase space, the 
z> = 3 fluctuations dominate the specific heat coeffi- 
cient 7. As a consequence, 7 changes its behavior at 
the (I)/ (II) crossover, £ ~ The thermal expansion, 
a, is at first sight also dominated by the z > = 3 mode. 
However, the repeated scattering of z < = 2 quantum fluc- 
tuations leads to a logarithmic dependence of the ther- 
mal expansion vertex. This in turn results in a sensitiv- 
ity of the thermal expansion on the (II) /(III) crossover 
too, and a therefore changes its critical behavior at both 
crossovers. Finally, the asymptotic critical behavior of 
the compressibility, k, is dominated by the z < = 2 mode 
and it is thus sensitive to the (II)/ (III) crossover only. 
To summarize, the critical thermodynamics knows about 
both crossovers, (I)/ (II) and (II)/ (III) in Fig. [5] and can- 
not be interpreted in terms of only a single dynamical 
exponent. 
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B. Outline of the paper 

In Sec. [TU we review the derivation of the effective 
Ginzburg-Landau theory of the Pomeranchuk instability 
in an isotropic metal of Ref. [TEl . In addition, we dis- 
cuss corrections to the shear susceptibility beyond the 
random-phase approximation (RPA), and we point out 
the importance of vertex corrections in order to main- 
tain the z < = 2 dynamics of the undamped shear mode. 
Readers mostly interested in the main results might di- 
rectly proceed to Sec. where the effective bosonic 
field theory is analyzed, and the critical thermodynamics 
is discussed. We end with a summary of the results in 
Sec. |IV] Detailed calculations arc relegated to the ap- 
pendices. 



II. POMERANCHUK INSTABILITY OF THE 
FERMI LIQUID 

We shortly review the derivation of the effective Ginz- 
burg-Landau theory for the Pomeranchuk instability of 
an isotropic metal by closely following Oganesyan et al<^ 
We limit ourselves to considerations in spatial dimensions 
d = 2. The starting point is a model of spinless fcrmions 
interacting via a quadrupolar interaction, 



S 



drd a x 



/t/ 



(3) 



where g 1 (k,iui n ) = —iuj n + £fc — /i is the fermionic 
Green's function. The operator Qij is defined by its 
Fourier transform with respect to the spatial coordinate, 

*£Qi3*x' -* *kQy[(k + k')/2]*k' where the quadrupo- 
lar momentum tensor is 



(4) 



with d = 2, kj = kj/fc and i,j £ {1,2}. The expecta- 
tion value Eij = (^'Qij^) is the quadrupole density and 
can be interpreted as the traceless part of a strain tensor 
representing the elastic shear modes of the Fermi sphere. 
With the help of the Hubbard-Stratonovich transforma- 
tion defined by the action 



S 



drd d x 



2F : 



(j^Ui-j - a i:j ( ^Q, A> 



(5) 



we can introduce the quadrupolar stress tensor that 
is real, symmetric, and traceless. Upon integrating out 
the fermionic degrees of freedom one obtains an effective 
elasticity theory for the Fermi liquid in terms of . 

General symmetry considerations require that its static 
Ginzburg-Landau energy functional has the form 



a Do 

12 



12 



(6) 



with the quadrupolar gradient tensor D^j =ddidj—SijV 2 . 
Note that generally in d = 2 a trace over a product of 
three quadrupolar tensors vanishes. In particular, a term 
tr{cr Da} is absent so that there is only a single Frank 
constant K on the quadratic level of the theory in the 
lowest-order gradient expansion. In the following, we will 
neglect the term with constant Ki as it is irrelevant for 
our considerations. We note in passing that the analyt- 
icity of the gradient expansion employed in Eq. ([6]) is ro- 
bust against interaction effects; it was argued in Ref. [12 
that the nonanalytic momentum dependences that are 
generated close to a ferromagnetic quantum critical point 
are absent for the Pomeranchuk instability. 

The quadrupolar gradient tensor also allows for a cou- 
pling of the shear modes to the compressional mode, i.e., 
the isotropic density fluctuation, 5n, in the Fermi liquid 
of the form Sntr{Da}. If this coupling is strong, it hy- 
bridizes the two different modes giving rise to low-energy 
excitations at some finite momentum. We will assume 
that this coupling is sufficiently weak such that we can 
neglect the compressional fluctuations in the following. 

The dynamics of the shear modes a t j are derived from 
an explicit calculation of the quadrupolar polarization of 
the Fermi liquid^ For this purpose, it is convenient to 
expand the shear field at a given momentum q and Mat- 
subara frequency Cl n , o"ij(q, if2„), into components lon- 
gitudinal and transversal to the quadrupolar momentum 
tensor, 

(Ty(q,«n n ) = 0|| q ,^„K^(q) + </>±q,in n -^(q). (7) 

The matrices E"' 1 - form a basis set for 2x2 real, sym- 
metric and traceless tensors and are defined with respect 
to the momentum direction q = q/q, 

E l (4) = -j= - Sij) = -j= {q % q 3 - piPj) , (8a) 



1 



Ejjiq) = (qtPj+Mj), 



(8b) 



where the unit vectors p and q are orthogonal with 
det(p, q) = 1. The basis set is normalized Ef^E^ = S a p, 
with a, j3 £ {||, _L}. The matrices of different basis sets 
transform according to 



£g(k) = L^(q) 

with the transformation matrix 

cos 2(f> sin 2(f> 
_ I - sin 2<p cos 2(f) 



(9) 



(10) 



where qk = cos <fi and pk = sin <j>. U describes a rota- 
tion but with angle 2<fi reflecting the invariance of the 
quadrupolar field under <j) — 7r rotations. 



A. Quadrupolar polarization of free fermions 

The quadrupolar polarization of the fermions, see 
Fig. [3£a), gives rise to a dynamical component of the 
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(a) 



(6) 





FIG. 3: Contributions to the quadrupolar polarization of the 
Fermi liquid. Solid lines are fermion propagators, the dashed 
lines are susceptibilities of shear fluctuations, and the vertex is 
given by the quadrupolar gradient tensor Qij . (a) polarization 
of free fermions, (b) and (c) lowest-order corrections. 



shear susceptibility tensor Xijki — (ffijCti), 

XoijkM^rv) = -j-^jki - n?. M (q,iQ„), (11) 

where t i3k i = \ (Sik$ji + 5ab 3 k - %4i) is the projection 
onto the 2d quadrupolar subspace and II is the lowest- 
order quadrupolar polarization of fermions, 



n°. H ( q ,ffi„) = - i QijMQkiOt) 



(12) 



In the basis ([8]) for the shear modes, the susceptibility 

X~ a l = E^XiluE^ reads 



X a p(q,in n ) = -— 6 a/3 -n%(q,iO n ) 



(13) 



where a, (3 E {||,_L}. The polarization at zero external 
frequency and momentum is given by the density of states 
11° ^ (0,0) = v8 a p. The criterion x~^(0,0) = for the 
Pomcranchuk instability to occur translates to^ 



(14) 



The dynamic part of the polarization is also diagonal in 
the basis ((8]). Its longitudinal part reduces to 



n|jii(q,ift„) 



f 27T dp 2 cos 2 2ip 
= —i\l n v \ — 

dyn Jq 2n iil n — VFq COS If 

w -2i/l^l. (15) 
v F q 



For the transversal part one obtains instead 



ni ± (q,in n ) 



dyn 



2n dip 2 sin 2 2ip 



-4, (ShL) 

\VFqJ 



o 27r ifl n — VFq cos ip 

2 



(16) 



The integral over the angle ip takes into account the pos- 
sible relative orientations of the mean Fermi velocity of 
the particle-hole pair and the momentum of the shear 
fluctuations, q; the p dependence of the weight distin- 
guishes the longitudinal from the transversal mode. The 





FIG. 4: Visualization of a quadrupolar fluctuation (black 
solid line) of the 2d isotropic Fermi sphere (gray solid line), 
(a) The longitudinal mode is damped by exciting particle- 
hole pairs close to the Fermi surface, Eq. ()15|1 . (b) For the 
transversal mode the momentum q of the fluctuation is close 
to a node such that there is not sufficient phase space for 
Landau damping, Eq. (|16p . 



approximations in Eqs. (jl5|) and (|16|) arc valid to lowest 
order in £l n /(vpq). 

The origin of the different dynamics of the two modes is 
illustrated in Fig. [U The quadrupolar fluctuation of the 
Fermi sphere cuts the undistorted sphere in four nodes. 
The nature of the Pomcranchuk dynamics depends on the 
relation between the location of these nodes and the di- 
rection of bosonic momentum q. In the low-energy limit 
io <C VFq, the bosonic momentum q is approximately 
tangential to the Fermi surface. Near the antinodes, see 
Fig. SJa), Landau damping is active and the longitudi- 
nal quadrupolar mode is damped by exciting particle- 
hole pairs, Eq. (fT5"]) . Near the nodes, on the other hand, 
the lack of phase space for exciting particle-hole pairs re- 
sults in an undamped dynamics for the transversal mode, 
Eq. IHl), see Fig.gfb). 

In the lowest order in frequency and momentum, the 
shear susceptibility [Eq. (fTB"]) ] is diagonal in the basis (jSJ), 



Xoa/3 



SapXoaai with the two components, 



Xo|n| fa' = r o + K q 2 + V 



Xo l ± (<l,in n ) =r a + Kq 2 + t 1 '— (17b) 



JOnl 

v F q ' 

n 2 

{v F qf 



(17a) 



where tq = — 1/(^2) — v, and 77 = ?//2 = 2v. Note that 
the Franck constant K, see Eq. ([6]), is not determined 
by the low-energy properties of the model [Eq. |2])] but 
depends on the details of the electron band structure. 

As advertised, from Eqs. (fl~7|) it follows that the dy- 
namics of the shear fluctuations are characterized by two 
different dynamical exponents that define the dispersion 
of the two eigenmodes, w ~ q z , at criticality r$ ps 0. 
Whereas the damped longitudinal mode has z> = 3, the 
dynamical exponent of the undamped transversal mode 
is instead z< — 2. Note that symmetry ensures that both 
modes have the same tuning parameter ro even after in- 
cluding interaction corrections and they become simul- 
taneously critical. The Pomcranchuk quantum critical 
point is thus a natural candidate to study the interplay 
of critical modes with different dynamics. 
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Corrections to the quadrupolar polarization of 
free fermions 



One may wonder whether the dynamics of the shear 
susceptibility [Eq. (|17|) ] is robust or gets modified in 
higher order loop corrections. The damped longitudinal 
shear fluctuations couple back to the fermions resulting 
in singular self-energy corrections^ The resulting non- 
Fermi liquid might screen the shear fluctuations in a dif- 
ferent manner than the Fermi liquid maybe leading to a 
modification of the effective shear fluctuation dynamics. 
It was speculated in Ref. that this feedback might al- 
ter the dynamics of the transversal mode and change its 
dynamical exponent. 

A detailed analysis of the ncxt-to-leading loop dia- 
grams in Figs.^b) and[3fc) is presented in Appendix |A"1 
We indeed find that by considering only the quadrupo- 
lar polarization diagram (b) containing the singular self- 
energy correction one would erroneously conclude that 
the dynamics of the transversal fluctuations gets modi- 
fied from z< = 2 to z< = 12/5, sec Eq. (|A12I) . How- 
ever, it turns out that the contribution from the singular 
self-energy in diagram (b) is just canceled by the vertex 
correction displayed in Fig. EJc). Similar cancelations of 
singular self-energy and vertex contributions in the elec- 
tron polarization are well known, for example, from the 
problem of electrons coupled to a Landau-damped gauge 
field. 33 i 34 ' 35 The root of this cancellation is the Ward 
identity deriving from particle number conservation.— 
Upon including both diagrams (b) and (c), we finally ar- 
rive at the result that the dynamics of the shear fluctua- 
tions [Eq. (jTTJ)] is robust with respect to the lowest-order 
corrections to the free fermion quadrupolar polarization. 
In the present problem of the Pomcranchuk instability, 
self-energy and vertex corrections thus have to be treated 
on the same footing in order to maintain the dynamics 
of the shear modes. Note that for this reason a simple 
version of Eliashberg theory, which neglects the vertex 
diagram Fig.[3]Jc), is not applicable. 

In the following, we will take the form of the dynamics 
for the shear fluctuations for granted that derive from 
the free- fermion polarization, Eqs. (jTTJ) . Using this dy- 
namics, we then consider the effective Ginzburg-Landau 
theory for the shear modes in the spirit of Hertz 4 



(a) 



(b) 



FIG. 5: One-loop corrections for the effective theory [Eq. (Tf8 
of the shear modes. 



term tr{er 3 } vanishes in d ~ 2. Expanding the shear 
modes in the basis ([5]) with respect to a fixed momentum 
direction, say E a (e y ) with det(ej,,e x ) = 1, the theory 
takes the form of a $ 4 theory of a two-component field 



sm = J 



drd d a 



Uv 1 

2 



4! 



(19) 



For the susceptibility Xi we use the free fermion form 
[Eq. (fT7)) ] derived in Sec. Ill Al It is a 2 x 2 matrix given 
by 



^ B qdiag{X|ML q ,in„-X±±q,in„}^ 



(20) 



where the transformation matrix Us „ is defined in 
Eq. (pTJj) . The longitudinal and transversal components 
of the susceptibility are 



x - | 1 (q, 2 a i )=r + 9 2 + ^>— , (21) 

(22) 



'1 

Q' 2 

— X I r> \ 2,2 "n 



with a conveniently rescaled momentum and tempera- 
ture; 7?> and ?7< are dimensionless (positive) parameters. 
Note that we could have chosen units such that cither 
V> 01 7 1< is unity so that one of them is in fact redun- 
dant] we deliberately keep here both explicitly for later 
convenience. After rescaling, the dimensionless length 
and temperature are measured in units of Kj v\ and 
[vpy/i/fK], respectively, where K is the Franck constant, 
v is the fermionic density of states, and vp is the Fermi 
velocity. The parameter r controls the distance to the 
Pomeranchuk transition and the interaction amplitude 



III. GINZBURG-LANDAU THEORY FOR 
POMERANCHUK QUANTUM CRITIC ALITY 

From the static energy functional [Eq. ([6])] and the 
shear susceptibility [Eq. (jTTJ)] we can construct an effec- 
tive Ginzburg-Landau theory for the shear modes, o^-, of 
the Fermi sphere, 



5 H -/ 



drd d 2 



(18) 



We limit ourself to spatial dimensions d = 2. As al- 
ready mentioned in the context of Eq. ([5]), a third-order 



A. Perturbative renormalization group at T = 

The transversal component of the shear modes with 
susceptibility xJ-L has dynamical exponent z < = 2 and 
thus an effective dimension d + z < equal to the upper 
critical dimension d+ = 4 of the <E> 4 theory [Eq. (fT9|) ]. 
Perturbation theory in the interaction u at temperature 
T = is therefore accompanied by logarithmic singular- 
ities which we sum up with the help of an one-loop RG 
treatment. 

Consider the perturbative corrections in Fig. [5] For 
their evaluations the angular averages of shear mode 
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propagators given in the Appendix [5] are needed. The 
one-loop self-energy diagram (a) leads to a shift of the 
tuning parameter r — > r + Sr, 



Sr 



u 1 

3/3 



(Xj|||q,lQ„ + X-L-Lq,in„) 



(23) 



q,O n 



Similarly, the vertex correction (b), u — > u + <5u, is given 
by 



<5u 



IT?; E [ 9 (4llq,<o» + xix q ,iQ„) (24) 

q,fin 



-2x||||q,in„X±±q,-tn„] 



As usual, the renormalization of the control parameter 
[Eq. (|23p] yields a nonuniversal shift, which depends on 
the UV cutoff of the integrals. In the following, we as- 
sume that such finite renormalizations have been already 
absorbed by appropriate counterterms.— In addition, the 
diagram (a) gives a universal albeit logarithmically diver- 
gent correction proportional to the control parameter, r, 
itself arising from the transversal z< = 2 quadrupolar 
mode, 



dSr 
dr 



u 1 

"3/3 



E Xj--Lq.if2„ 



■ A 

247rr,<' ° g I 



(25) 



In the last expression, we evaluated the integrals at T = 
with a sharp momentum cutoff A and kept only the lead- 
ing logarithmic term. Similarly, we get for the vertex 
correction 



Su 



2v* 



E ^J--Lq,in„ 



3u 2 



q,n„ 



327TT7< 



■log 



A 



(26) 



Equations. ([23)1 and (|2"6"|) are the only logarithmic sin- 
gularities encountered in lowest-order perturbation the- 
ory. In order to treat these singularities, we apply the 
renormalization group and perturbativcly integrate out 
modes within a momentum shell (A/6, A) with log 6 >C 1. 
After such a single RG step, the momentum cutoff of 
the perturbatively corrected theory is restored to A by 
rescaling momenta q — ► q/b and frequencies ui — > uo/b z 
with some dynamical exponent z. As we treat a theory 
with multiple dynamical scales, one may wonder what 
kind of value to choose for z in this RG process. It is in- 
structive to leave the value for z here unspecified in order 
to discuss later possible choices. Of course, the result for 
physical observables will turn out to be independent of 
the precise choice for z. Using the renormalization-group 
conditions that the prefactor in front of the momentum 
dependence, q 2 , of the susceptibilities, Eqs. (|2"Tj) and (|22|) . 
remains unity, the pcrturbative rcnormalizcd parameters 
can then be read off. The scaling of the dynamical pa- 
rameters, rj < and ry>, at one-loop order is simply obtained 



from power counting 



d log b 
d log b 



(z> - z)?7>, 
(z< - z)t7<, 



(27a) 
(27b) 



where again z< = 2 and z> = 3. The one-loop cor- 
rections to tree-level scaling for the control parameter, 
r, and the interaction, u, at T = follow directly from 
Eqs. (J2SJ) and (J2BJ, 



dr 
d log b 

du 
d log b 



(4- d — z)u- 



327r?7< 



(28a) 
(28b) 



with d = 2. These RG equations are characteristic for 
the universality class of the Pomeranchuk transition and 
are an important result of this work. 

Let us remark on the ambiguity with respect to the 
choice of z, that determines the tree-level scaling in 
Eqs. ((17|) and (25]). First of all, note that for all choices 
of z > 2 the tree-level scaling suggests that the quar- 
tic coupling u as well as the dynamical parameter r] < are 
irrelevant Nevertheless, we know from the perturbative 
results, Eqs. (f2"S"|) and (|2"6"|) , that there exists a nontrivial 
marginal operator. Indeed, closer inspection of the RG 
equations readily shows that the quartic coupling enters 
the RG flow only in the combination 



u 



Correspondingly, the Eqs. 
Eq. (j27b|) to give 

dr 



(29) 



can be combined with 



d log b 

du 
d log b 



2 u I r, 

3 ~2 
u . 

32tt 



(30a) 
(30b) 



The seemingly irrelevant parameters u and rj < thus con- 
spire to yield in fact a marginal operator with amplitude 
u = u/rj < . Note, in particular, that Eqs. (|30|) are inde- 
pendent of the specific choice of z. 

The freedom to choose the exponent z can now be ex- 
ploited to simplify the RG flow in Eqs. (|27|) . There are 
two evident choices, z — z < or z = z > , that eliminates 
the flow of one redundant scaling variable, i.e., cither 
r) < or n > , respectively. For z = z < the remaining flow 
of the dynamical parameter n > is unbounded, n > — > oo 
whereas for z — z> the other parameter flows to zero, 
rj < — * 0. That means that depending on the choice of z 
either l/r/> or n < is irrelevant in the RG sense. Taken 
together, the RG flow has either the critical fixed point 
r = u = l/?7> = or r = u = ij < =0, respectively. 
However, note that both fixed points are only of limited 
use when computing physical properties. The reason is 
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that both, l/?7> and 77<, are in fact dangerously irrelevant 
for the respective choices of z. When considering, e.g., 
the critical thermodynamics, it turns out that neither 
l/?7> — > nor r\ < — > is a well-defined limit, which is the 
very definition of a dangerously irrelevant variable^ see 



Appendix IC II Nevertheless, let us stress that any choice 
of z for the RG yields, of course, the same result when 
computing critical properties of physical observablcs. 

The Eqs. (|30[) arc readily solved. The running coupling 
u(b) is given by 



u(b) = 



32tt 



1 



3 log(A6/A) 



(31) 



with A = Ae 327r? '</ (3tl) . The scale dependence of the 
coupling u = u/rj < can be interpreted as an effective 
quartic vertex, 



u e $(q,u,T,£) 
32irr/ < 



(32) 



log(A/max{g, W V*< , T l / Z < , !/£}) : 



that depends on momentum, frequency, temperature and 
correlation length £. The full quartic vertex generally 
depends on four momenta and frequencies. However, to 
logarithmic accuracy we can take q = max-j^} and 
lu = maxjw,}. It is important to note that it is the 
dynamical exponent of the transversal mode z < = 2 ac- 
cording to which the momentum scale is translated to 
energies in Eq. (f3"2"| as it is this very same mode that 
induces the logarithmic dependence of the vertex. 



This logarithmic enhancement of the correlation length 
was already advertised in Eq. (JT|). The exponent 4/9 is 
characteristic for the Pomeranchuk universality class. It 
is instructive to compare this exponent with the value ob- 
tained from a 0(A^)-symmetric $ theory with a propaga- 
tor containing only the z< = 2 dynamics, Xa/3 = ^a0X±j_ 
and a, j3 = 1, . . . , N. In such a case, one would obtain 
the exponent (N + 2)/(N + 8) instead of 4/9. In par- 
ticular, note that the Ising (N = 1) value 3/9 as well as 
the XY (N = 2) exponent 4/10 differ from the value we 
obtain for the Pomeranchuk transition. The difference to 
the Ising exponent 3/9 clearly indicates that the model 
[Eq. (fT9|) ] cannot effectively be reduced to a $ 4 theory 
of a single-component field with dynamics z < = 2 as one 
might naively expect. The origin of the anomalous value 
4/9 of the exponent in Eq. (f3"3"| can be traced back to the 
result of angular averages over shear- mode propagators, 
Eqs. (|B1[) . and thus reflects the topology of the shear 
field. 

At finite temperatures, the correlation length also ac- 
quires a T dependence. In the low-temperature regime (I) 
of the phase diagram in Fig. [2J this T dependence is only 
a sub-leading correction to Eq. (f33)) and can be computed 
within renormalized perturbation theory from the dia- 
gram in Fig. 03(a). It is the thermal occupation of z> =3 
fluctuations that yields the dominant T correction, which 
is of Fermi-liquid type, £ -2 — S,t=o ~ @( u cttT 2 {; 3 ), where 
Moff ~ m(A£) is the quartic coupling [Eq. (|3"Tj) ] at the scale 
b = A£. 



B. Correlation length 

We compute the dependence of the correlation length 
£ on the control parameter, r, and temperature, T. We 
distinguish between the regimes (I), (II), and (III) of the 
phase diagram in Fig. O Whereas in regimes (I) and 
(II) the correlation length is determined by the zero- 
temperature flow of the control parameter, Eq. (|28a|l . 
it is dominated by temperature, T, in regime (III). Re- 
markably, this latter T dependence is found to be univer- 
sal independent of the quartic interaction u. This uni- 
versality arises from a peculiar interplay of logarithmic 
singularities of different origin. 



1. Correlation length in regime (I) 

Consider first the correlation length at temperature 
T = 0, where it is determined by the RG flow of the con- 
trol parameter [Eq. (|28a|) ]. With the help of expression 
(f3Tj) for the running quartic coupling, we can solve for the 
running r(b) from which follows the correlation length at 



zero temperature £ \t= 



r{b)/b 2 



IT=0 



4/9 ' 



(33) 



2. Correlation length in regimes (II) and (III) 

The temperature dependence of the correlation length 
in regimes (II) and (III) arises from a subtle interplay 
of logarithmic IR singularities driven by quantum and 
classical fluctuations. In order to illustrate this peculiar 
low-energy structure of the theory, consider first the T 
dependence of the control parameter deriving from the 
perturbative one-loop correction (a) in Fig. [5l 



5r\ 



T=0 



dqq 
2tt 



dui / oj 
— coth 1 

7T V 2T 



x Ii: n{x!!!l+io+xii+io} 



(34) 



As before, it turns out that it is the contribution of the 
longitudinal susceptibility x"" that gives the leading T 
dependence in regimes (II) and (III). Furthermore, in 
d = 2 the integral [Eq. (|34j) ] is dominated by the limit 
of small frequencies co in these regimes, i.e., the classical 
limit of the hyperbolic coth function corresponding to 
the contribution of the Matsubara zero mode. Generally, 
the classical limit is applicable whenever the frequency 
is smaller than temperature, lo < T. Using the scaling 
relation, ui ~ q z > this translates to a regime of momenta, 



1 < Ct i w ith the thermal length £^ ~ T 1 / z >, where 
the associated fluctuations have a classical character, see 
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Fig. m 



Sr - Sr\ T=0 



uT 
~3~ 



dqq 
2tt 



X llll (q,^ = 0). (35) 



In the classical regime, the renormalized susceptibil- 
ity reduces to Xyii 1 ^) = ( 2 + <Z 2 ' an d the momentum 
integral in Eq. (|35[) will be cut off by the inverse cor- 
relation length l/£ at small q. Anticipating that the 
correlation length itself is cither given by £ ~ 
in regime (III) or approximately by Eq. (|33p in regime 
(II). the momentum integral in Eq. (|35[) extends over 
a parametrically large range l/£ < q < T 1 / 3 , where it 
is logarithmically enhanced. However, in this momen- 
tum range the transversal mode still possesses its quan- 
tum character, see Fig. Q] As we showed in the previous 
section, for momenta larger than the inverse correlation 
length, q > l/£, the vertex acquires an effective scale 
dependence due to the repeated scattering of quantum 
fluctuations in the transversal channel, see Eq. ([32]). As 
the integral in Eq. (|35|) is logarithmically enhanced, the 
temperature correction turns out to be sensitive to this 
scale dependence of the vertex which becomes apparent 
in higher-order loop corrections. The scale dependent 
vertex [Eq. ([32]) ] is thus probed by the classical fluc- 
tuations in the longitudinal channel with momenta be- 
longing to the extended quantum-to-classical crossover 
regime (T 1 / 2 , T 1 / 3 ) of Fig.CD In order to determine the 
correlation length it is therefore insufficient to limit one- 
self to the lowest-order diagram (a) . 

Instead, we apply the RG generalized to finite T, see 
Appendix [C] When the cutoff is lowered down to a value 
on the order of the inverse thermal length £^ 1 corre- 
sponding to a scale b ~ A£^, the contribution arising 
from the integral [Eq. (f3"S")) ] translates to a modification 
of the zcro-tempcraturc flow in Eq. (|28a|) . The classical 
longitudinal fluctuations give rise to an accelerated flow 
of the running control parameter r(b) setting in at the 
scale b ~ , 



_i_^ r + i_ uT e 6 _ Af > (36) 



d log b 



where Q x — 1 if x > and otherwise. The flow of 
the quartic coupling is still given by Eq. (I28b[) and the 
temperature has the scale dependence T(b) = Tb z . Solv- 
ing for r(b) and identifying the correlation length with 
£~ 2 = r(b)/b 2 \b = \£, we get the following implicit equa- 
tion for £ valid in regimes (II) and (III) of the phase 
diagram in Fig. [21 



/ log A/ A A 
V logA£ J 



4/9 



A V< T 1 - 



logA£> 
logA£ 



4/9N 



(37) 



with £^ ~ T" 1 / 3 . The first term dominates in regime 
(II) and yields a value for £ that coincides with expres- 
sion ([3"3"| . The second term originates from the accel- 
erated RG flow due to the classical fluctuations in the 



longitudinal channel. This acceleration acts on scales 
in the window < b < A£, which can be identified 
with the momentum range in Fig. ([T]) where the quan- 
tum regime of the transversal mode overlaps with the 
classical regime of the longitudinal fluctuations. This T 
dependence dominates the correlation length in the quan- 
tum critical regime (III). 

At criticality r = 0, Eq. ([3~7| is easily solved itcrativcly 
and we get a universal asymptotic behavior at low tem- 
peratures, 



r 



lr=0 



Cr?<T, C 



4/9' 



(38) 



This reduces to Eq. (J2J) in Sec. II Al after exploiting the 
freedom to choose temperature units such that rj < = I, 
cf. the remark below Eq. ([2"2"| . Note that the temper- 
ature dependence of £ is here independent of the bare 
value of the quartic coupling u. It enters the implicit 
Eq. ([3T)) via A defined below Eq. ([3"Tj) but drops out in 
the asymptotic limit T — > 0. This universal tempera- 
ture dependence originates from an intricate interplay of 
quantum fluctuations in the transversal channel and the 
classical shear fluctuations in the longitudinal channel. 



3. Ginzburg temperature 

Our perturbative treatment finally breaks down when 
classical loop corrections start to dominate. This is the 
case when the Ginzburg parameter is on the order of one, 
Q = u c ffT/£~ 2 , where u e s ~ w(A£) is the running vertex 
[Eq. (|3T])] at the scale b = A£. The criterion Q ~ 1 iden- 
tifies the Ginzburg temperature Tq(t). In the asymp- 
totic low-temperature limit, the Ginzburg temperature 
is given by 



T G (r) 



and is shown in Fig. [2j 



log 



4/9 



(39) 



C. Critical thermodynamics 

The presence of two dynamical scales gives rise to two 
crossover lines in the phase diagram, see Fig. [2l In the 
present section, we investigate how these crossovers affect 
the critical thermodynamics. 

The critical free energy, J- CI , depends on the control pa- 
rameter, r, and temperature, T. We consider the thermo- 
dynamic quantities that correspond to the three second- 
order derivatives of J-" cr with respect to r and T. The 
derivative, 7 cr = —d 2 J 7 cl /dT 2 is the specific heat coeffi- 
cient. We discuss the derivatives with respect to r in the 
language of a pressure-tuned quantum phase transition. 
If the control parameter depends smoothly on pressure, 
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r = r(p), derivatives with respect to r effectively corre- 
spond to derivatives with respect to Consequently, 
the quantity a cr = d 2 J- CI /(dTdr) can be identified with 
the thermal expansion and K a = —d 2 T cr /dr 2 parallels 
the critical contribution to the compressibility. 

It is instructive to discuss the asymptotic behavior 
close to criticality within a renormalized Gaussian model, 
where the free energy is given by 

d 2 q 1 



1 



ren. Gauss 



(2n) 2 (3 ^ 



logdetx (q,ifl n )\r 



(40) 

with the susceptibility tensor of Eq. (|2T))) but with the 
control parameter replaced by the correlation length, 
r — > t 2 - The approximate expression (|40"|) for the free 
energy accounts for most of the critical behavior except 
the thermal expansion in regimes (II) and (III), and the 
compressibility in regime (I) of Fig. O In order to ob- 
tain the correct asymptotic behavior of the latter, we 
have to resort to a proper RG improved treatment. The 
quantum fluctuations of the z < = 2 mode lead here to 
scale-dependent vertices that are beyond the renormal- 
ized Gaussian form of Eq. (|40|) . Below, we only present 
the result of the asymptotic critical behavior and refer 
the reader to Appendix [C] for its derivation. 



1. Specific heat coefficient 

The critical specific heat, 7 cr = —d 2 J 7 cr /dT 2 , is domi- 
nated by the damped z > = 3 shear fluctuations because 
their available phase space is larger. Their contribution 
to the specific heat can be derived within the renormal- 
ized Gaussian model with a free energy given by Eq. (|40[) . 
In the low-temperature regime (I) of Fig. [2j the criti- 
cal specific heat is temperature independent in the limit 
T ^ 0, 



7cr ~ £ 



r l/2 



2/9 



CD- 



(41) 



In the overlap (II) and quantum critical regime (III), we 
get 

7cr ~T- 1/3 (II) and (III). (42) 

The undamped z < = 2 mode adds a sub-leading correc- 
tion, which is on the order 0(T£ 2 ) in regimes (I) and 
(II), and 0(1) in regime (III). Only the sublcading cor- 
rections to 7 cr distinguish the regimes (II) and (III) . Note 
that the z< = 2 mode, nevertheless, leaves its traces in 
7 cr via the logarithmic rcnormalization entering the cor- 
relation length in regime (I). 



2. Thermal expansion 

The dependences of the thermal expansion, defined as 
a cr = dJ- CI / (drdT), is more delicate. It is again the large 



phase space of the z > = 3 mode that is the main driving 
force. In the low-T regime (I), a CI can still be understood 
within the renormalized Gaussian model [Eq. (140]) ]. 



3^ 
dr 



T 

^3/2 



log- 



A 



2/!) 



(I). (43) 



This leading contribution is attributed to the z > = 3 
part of Eq. (|40p . The z < = 2 mode again enters here 
indirectly in the form of the logarithmically renormalized 
correlation length. 

Interestingly, the critical thermal expansion, a cx , with- 
in the regimes (II) and (III) is more involved. Let us 
first discuss the result one would obtain from the renor- 
malized Gaussian model [Eq. (|40|) ]. In d = 2, it is the 
classical Matsubara zero mode that gives the dominating 
contribution to a C r, 



dq 



at 



t 



di 



(44) 



We only retained the part due to the classical z > = 3 
mode here; the additive classical z < = 2 contribution 
is cut off at the smaller momentum scale £^ 1 instead 
and turns out to be subleading. The derivative df^ 2 /dr 
can be interpreted as an effective vertex that is, in gen- 
eral, scale dependent due to the RG flow of the control 
parameter discussed in the context of Eq. (|28a[) . As 
the momentum integral in Eq. ([44]) covers the overlap 
regime (l/£, l/£^), where the z < = 2 fluctuations arc 
still quantum, cf. Fig. Q] this scale dependence is fully 
developed in the important integration range and has to 
be taken into account. Similarly as for the correlation 
length, the quantum fluctuations of the z< = 2 mode 
here give rise to a scale dependent vertex that is probed 
by the classical fluctuations of the z > = 3 mode. The 
critical thermal expansion is thus beyond renormalized 
Gaussian but can be obtained from Eq. (|4"4"|) by replacing 
d^ 2 /dr — > l/(log(A/g)) 4 / 9 . So we finally get to loga- 
rithmic accuracy 



dq 



1 



q [log(A/g)]4/9 



(II) and (III) (45) 



whose asymptotic limit depends on whether the correla- 
tion length is temperature dominated, (III), or r domi- 
nated, (II). To summarize, we have 



Or 



T 

r 3/2 

(log 



( log 7?) 



2/9 



5/9 



(D 
(ii) 
(in) 



(46) 



The thermal expansion thus changes its behavior at both 
crossovers in Fig. O 
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3. Compressibility 

The scale-dependent vertex in Eq. ([45]) is also impor- 
tant for the critical compressibility, K cr = —d 2 J 7 cr /dr 2 . 
In the low-temperature (I) and overlap regime (II), the 
Gaussian free energy [Eq. (|40[) ] predicts a logarithmic 
divergence of the critical compressibility due to the un- 
damped 2< = 2 shear fluctuations 

The correct behavior again follows upon RG improvement 
by replacing the vertex d£~ 2 /dr — ► l/[log(A/fc)] 4 / 9 , 

( A\ 1/9 
K cr ~ (log-^=) (I) and (II). (48) 

Alternatively, this result can be derived directly from the 
non-trivial RG flow of the critical free energy as explained 
in Appendix [C] Note that the logarithmic divergence of 
the compressibility as r — ► + indicates that a coupling 
of other degrees of freedom to the square of the shear 
held tr{er 2 } might be nonperturbative^ 

In the quantum critical regime (III), on the other 
hand, the correct asymptotics can again be obtained from 
Eq. (001), 

/ A V 8/9 
K cr ~(log-^J (III). (49) 

Both modes contribute to this logarithmic dependence 
of the compressibility. While the z > = 3 fluctuations 
contribute only with their classical component, i.e., Mat- 
subara zero mode, remarkably, all Matsubara frequencies 
of the z < — 2 fluctuations add to the logarithmic depen- 
dence [Eq. (j4*9")) ] because £ ~ £^ in regime (III). 

So we find that the asymptotic behavior of the critical 
compressibility is only sensitive to the (II) /(III) cross- 
over. 

IV. SUMMARY 

We considered the problem of multiscale quantum 
criticality by studying the effective bosonic model of 
a quadrupolar Pomeranchuk instability in an isotropic 
metal in spatial dimension d = 2. This theory consists of 
two critical bosonic shear modes of the Fermi sphere char- 
acterized by different dynamics, resulting in a quantum 
criticality with multiple scales. While the longitudinal 
shear fluctuations are Landau damped and have a crit- 
ical dynamical exponent z > = 3, the transversal modes 
are undamped with z < = 2J£ 

First, we studied the case of zero temperature. There 
we found that the low-energy properties are governed by 
the mode with the smaller effective dimension d + £< , 
i.e., the transversal shear fluctuations. Since that effec- 
tive dimension is equal to the upper critical dimension, 



e?+ = 4, of the effective 4> 4 theory, its self-interactions 
give rise to logarithmic singularities in perturbation the- 
ory and, as a consequence, to a nontrivial RG flow of the 
theory, Eqs. (|28[) . The interaction amplitude is found to 
be marginally irrelevant in the RG sense so that the low- 
energy quantum critical fixed point is Gaussian. Never- 
theless, as usual the slow flow of the interaction ampli- 
tude toward weak coupling gives rise to qualitative cor- 
rections to Gaussian behavior. These corrections con- 
sist of, e.g., a logarithmic enhancement of the correlation 
length, Eq. (f3"3"| and a logarithmically diverging com- 
pressibility, Eq. (|4"8")) . whose functional form is charac- 
teristic for the Pomeranchuk universality class. 

Second, we investigated the critical properties at finite 
temperatures. A salient feature of quantum criticality, 
in general, is the quantum-to-classical crossoveri££i22 
associated with the thermal length, £t ~ T~ x l z . In 
the language of crossover theory^ the renormalization- 
group flow of a quantum critical model on short scales 
is first governed by the primary quantum critical, T = 0, 
fixed point, but at a scale on the order of £t the flow 
crosses over toward the secondary classical critical fixed 
point, that controls it for large scales. In the presence of 
multiple dynamical scales, we find instead a quantum-to- 
classical crossover regime, which extends over a paramet- 
rically large range bounded by the two thermal lengths, 
£> ~ T~ x l z > and ^ ~ T _1 / z <. In this extended cross- 
over regime, the transversal shear fluctuations (z< = 2) 
still possess a quantum-mechanical character and drive 
the theory toward the primary quantum fixed point 
whereas the longitudinal mode (z> = 3) is already clas- 
sical and competes to push the theory toward the sec- 
ondary classical fixed point, see also Fig. [TJ 

We find that this extended quantum-to-classical cross- 
over determines the critical properties in the quantum 
critical regime (III) of the phase diagram in Fig. O Re- 
markably, the correlation length, £, depends here on tem- 
perature in a universal manner, Eq. (|38p . This univer- 
sality emerges via a peculiar interplay of logarithmic sin- 
gularities of different origin. These singularities arise as 
the transversal quantum fluctuations have an effective di- 
mension d + z < — d^ equal to the upper, <i+ = 4, and 
the longitudinal classical fluctuations have an effective 
dimension d = d~ equal to the lower critical dimension, 
d~ = 2, of the respective $ 4 theories. We showed that 
the interaction of these modes results in a T-dependent 
boost in the RG flow of the control parameter, Eq. (|36[) . 
that finally determines the T dependence of £. 

The two thermal lengths, £^ and also identify 
crossover lines in the phase diagram Fig. [5] where the 
critical thermodynamics changes qualitatively. The sen- 
sitivity on these crossovers, however, differs for different 
thermodynamic quantities depending on whether they 
are dominated by the z > — 3 or the z < = 2 fluctuations. 
As a consequence, an interpretation in terms of scaling 
relations involving only a single dynamical exponent is, 
in general, not possible. 

To summarize, wc found a series of phenomena arising 
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from the presence of coexisting dynamics, which might be 
generic for quantum critical systems with multiple scales. 
So we identified (a) an extended quantum-to-classical 
crossover as summarized in Fig. [TJ with an intermediate 
momentum regime where quantum and classical fluctua- 
tions coexist and their interaction profoundly affects the 
critical behavior. The multiple scales are (b) reflected in 
the crossover lines in the phase diagram, see Fig. O We 
found that (c) certain type of thermodynamic quantities 
are governed by one or the other dynamics so that (d) 
simple scaling relations in terms of a single dynamical 
exponent are not applicable. Our study shows that mul- 
tiple scales at quantum criticality can yield a rich phase 
diagram, complex critical properties, and even emerging 
universality. 
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APPENDIX A: EFFECTIVE FERMION-BOSON 
THEORY FOR THE POMERANCHUK 
QUANTUM PHASE TRANSITION 

In order to study the shear susceptibility close to the 
Pomcranchuk quantum critical point beyond the free- 
fermion form, viz., RPA, see Eqs. (fT7|) . we consider the 
following fcrmion-boson model: 



(Al) 



drd d 



ijkl 



k + q 



k + q 1 + q 



k kl k-q « k 

(a) 



k + q' 



[b) 



FIG. 6: (a) Fermion self-energy and (b) vertex correction of 
the model [Eq. (|A1|| ]. For clarity, we dropped the frequency 
labeling. 



1. Fermion self-energy 

We calculate the fermionic self-energy, £ = g _1 — g$ . 
In second order in the coupling A, see Fig. EJa), it reads 



x XiiH(q,^n)fl , o(k-q,ia; Tl -jn n ). (A2) 

Taking into account that q is small compared to the 
Fermi wave vector this simplifies to 



£(k,iw n ) «-A 2 Q^(k)Q w (k) 



1 q,n„ 



(A3) 



-1 



iw n - iQ n 



Jfc + v F q 



with vp = w^k The important contribution to the self- 
energy originates from the longitudinal part of the sus- 
ceptibility 



£(k,^ n )«A 2 ±£ r^ XM ^iQ n ) (A4) 



2iT 



2 cos 2 2tp 

x / dip . 

Jo IUJ„ - lYL n - £fc + V F q COS If 

Close to the Fermi surface, this evaluates to 



where go is the fermionic Green's function given below 
Eq. ([2]), Xo is the shear susceptibility [Eq. (fT7|) ]. and the 



vertex is T 



ijkl 



lyfezA with the coupling constant A. 



The bosonic self-energy in second order in A, see 
Fig. E^a), will only shift the bare parameters of the sus- 
ceptibilty xo(q, i&n) but mantains the form of its U, n and 
q dependence. To indicate this renormalization, we will 
drop the index and consider the parameters of \ as 
having effective values that have to be determined self- 
consistcntly. The corrections in higher order C(A 4 ) to 
the quadrupolar polarization in Figs. [31(b) and [3tc) in- 
volve the lowest-order fermionic self-energy and vertex 
correction, respectively. We consider them first before 
turning to the discussion of the corrections to the shear 
susceptibility 



|w„| CO 

£(iw„) = -i-^— -signw„ J dfl J dqx\\\\(q,ity (A5) 

o o 

A 2 J ^ /3 i|a;„| 2 /3 s ig na , n if | Wn | » | r |3/2 

if \uj n \ < r 3 / 2 



2tt 



1/3 1 / 

where u> n = — -y= ^f- ■ Note that the leading-order 
contribution only involves the frequency dependence. 

The self-energy of the form (|A5[) is well-known from 
the problem of electrons coupled to a Landau-damped 
gauge fiel d 33 ' 34 ! 35 and ferromagnetic quantum criticality 
in metalsi^ 2 - 
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2. Vertex correction 



The lowest order vertex correction, ST 1 ^ 1 = yw^—Tq , 
shown in Fig. [S^b) reads 



50k+q',iw„+if2' n ffOk+q+q' ,io!„ +i£l„ +if2; 



(A6) 



where we used already that the bosonic momenta, q and 
q', are small compared to the fermionic momentum k. 
Applying the algebraic identity AB = — a AZb-i ^° * ne 
product of fermionic Green's functions, the expression for 
the vertex correction can be rewritten as 

«tt>>^n)« (A7) 

_ . Fi jk iCk + q, iw n + i£l n ) - F ijk i(\s., iu> n ) 
iQ n - v F q 

Here, we introduced the auxiliary function 
F ijkl {k,iu> n ) = -A 2 Q„-(k)Q mn (k) (A8) 



k+q', 

I'M 



Identifying Q mn (k) = V2-Emn(k) with the basis matrix 
defined in Eq. ([5]) and using the transformations [Eq. @] 
this simplifies to 

FijuOciUn) =-2A 2 4-(k)i ^ .9o k+q ',^ + ^ (A9) 

x (^wWcos^^xnilq'.in; 
+£^(k) sin 2 2cf>xxx< i >,in 

where kq' = cos</>. In the low-energy limit, this expres- 
sion reduces to Fijki (k, f=a 25y(k^i?|j(k)£(wt> n ) with 
the self-energy of Eq. (|A5j) . From Eq. (1A7I) thus follows 
the asymptotic relation for the vertex correction 



^„(q^)«-A4(k)4(k). 



J iuj n -\-i£l n 



itt n - v F q 



(A10) 



This relation can be identified as an asymptotic Ward 
identity^ deriving from fermion number conservation. 
The interaction vertex T^ kl in Eq. (|A1[) locally couples 
the fermion density ~^ k+q ^k to the projection of the shear 
field onto the fermionic quadrupolar momentum tensor, 

i.e., Elfaaijfatiln). (Note that Q y (k) « v^E^Ck).) 
This projection is thus a field conjugate to the density 
fluctuations. In other words, the projection of the vertex 
p*J fci _g;|^(k)£ , ^(k) is just a local density-density interac- 
tion and obeys the corresponding Ward identity, from 
which follows Eq. (|A10[) in the low-energy limits 



3. Interaction corrections to the quadrupolar 
polarization 

We now consider the modification to the bosonic prop- 
agator x hi fourth order in the fermion-boson interaction 
A due to the diagrams (b) and (c) in Fig.[3j In particular, 
we would like to check the stability of the dynamical part 
of x and its associated dynamical exponents. The strong 
longitudinal z> = 3 fluctuation arc responsible for the 
anomalous fermionic self-energy [Eq. (|A5[) ] of non-Fermi 
liquid form. This self-energy correction might feed back 
into the quadrupolar polarization of the electrons mod- 
ifying the dynamics of the perpendicular mode. So we 
focus in the following on the corrections to the perpen- 
dicular dynamics and show that its dynamical exponent 
z< = 2 is in fact robust against this feedback effect. 

First, consider the contribution of diagram (b) to the 
perpendicular mode, 



Aiv\ 2 
(VFq) 2 



duj (Uliuj + id, 



E(iw)) . 
(All) 



Using the expression (|A5|) for the self-energy, at critical- 
ity, r = 0, this reduces to 



<m£l(q,io,) 



24 AV i /3 10nl 5/3 

'5^ W ° (v F q) 2 ' 



(A12) 



The resulting dynamics seems to dominate over the low- 
est order result for n^j_ given in Eq. (fTTj)) . Taking into 
account only diagram (b), one might naively expect a new 
dynamical exponent for the perpendicular mode given by 
z< = 12/5 > 2. 

However, it turns out that the leading contribution of 
diagram (b) given in Eq. (j A12|) is exactly compensated 
by the diagram (c) that includes the vertex correction, 



(A13) 



x Qmn(k)Qfc;(k)gok+ q ,iw„+iQ„5ok,iw„ 

The cancelation of the leading contributions of diagrams 
(b) and (c) becomes apparent and is easily verified upon 
using the asymptotic identity for the vertex correction 
[Eq. {EE)]. 

A similar cancellation of singular sclf-cncrgy and ver- 
tex correction in the electron polarization due to the 
particle number conservation was discussed before in the 
context of the electron-gauge-field problem.— In particu- 
lar, note that a Eliashbcrg-type of theory would disregard 
the diagram (c) in Fig. [3] containing the vertex correc- 
tion and would erroneously conclude that the longitudi- 
nal z > = 3 mode couples back and modifies the dynamics 
of the perpendicular mode from z < = 2 to z < = 12/5. 
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APPENDIX B: ANGULAR AVERAGES OVER 
THE SHEAR MODE PROPAGATOR 

For the computations of the perturbative corrections 
in Sec. IIII Al the following angular averages are needed 



^ Xq% n = £ (Xllllq.ifin + X±-Lq,ifi„) 8 a p, (Bla) 



1 

2tt V 



(Bib) 



1 2 

g (X||||q,ifi„ - X_L_Lq,iO„) (^07^/35 + <5q5^ 7 ) 
+ g (x|||q, l r2„+xi± q ,ifi„+ 6 X||||q, J :0„X-L-Lq, t n„) (5 Q/3 < 57>5- 



APPENDIX C: RENORMALIZATION GROUP AT 
FINITE TEMPERATURES 

In this section, wc extend the RG treatment of the 
effective model [Eq. ([THjl ] to finite temperatures in order 
to treat the quantum-to-classical crossover. Following 
Ref. [H, we apply standard crossover theory^ to express 
the free energy as a line integral along an RG trajectory, 



T 



dXe -(d+z)x ^ f a {T(e x ),r(e x ), Va (e x )), 

«£{<,>} 

(CI) 



where d = 2. We follow the considerations in Sec. IIII Al 
and leave the exponent z unspecified in order to demon- 
strate the equivalence of different implementations of the 
RG procedure. The running temperature, T(b), control 
parameter, r(b), and dynamical parameters, rj a (b) are 
defined below. The integration kernels are given by 



fa 



A- 



8 1 



d 2 q 



OA 2 J (2tt) 2 /3 



(C2) 



where the identification _L = < and || = > for the index a 
applies, and the susceptibilities are defined in Eqs. I|2ip 
and (|22p . The momentum integral and the sum over 
bosonic Matsubara frequencies is understood to be reg- 
ularized with some UV cutoff A. Replacing the Matsub- 
ara sum with an integral by analytic continuation and 
employing a hard cutoff rcgularization, the function f a 
becomes 



f« = A 



d K 2 



dqq 



r 



dhJ , / u> \ 
— coth — 
2tt \2T) 



dA 2 Jo j_ A 

x Imlog ( X ~l(q,u + iO)) 



(C3) 



with K2 = l/(27r) and a principal- value integral over 
frequencies uj. Note that we choose to take the frequency 
cutoff to be A z . It is convenient to separate the functions 
f a into nonuniversal and universal parts, 

f a (T,r,r) a ) = fo, a (r,r) a ) + foo, a (T,r,r) a ). (C4) 



The nonuniversal parts, /o, Q , are given by 



fo<(r, n<) 
fo>(r, V>) 



= /<(0,0,7 ?< ) + /f 1 ' 0) (0,0,7 7< )r, 
= />(0,0,7 ?> ) + /( 04 ' 0) (0,0,7 7> )r 



(C5) 
(C6) 



1 



where f a °' n ^ is the nth derivatives with respect to the 
second argument. They contribute considerably only at 
the initial stage of the RG flow, i.e., for small value of A 
in Eq. (|C1[) , and give rise only to finite shifts of the bare 
values of the theory that can be absorbed by appropri- 
ate countcrterms. The universal parts of the integration 
kernels, /oo, Q , read explicitly (for A — > 00) 



/oo,<(T, r, ?7<) 



foc,>(T, r, r?>) 



1 



32ni] < 
A2 r oo 

-in 



(C7) 



du> 



coth ( — ) - 1 
\2TJ 



r 3 A 



157r 2 ?7>(r + A 2 ) 
A 2 



47T 2 Jo 

x arctan 



du) 
A(r + A 2 ) 



coth (—) - 1 
\2TJ 



(C8) 



The first terms on the right-hand side give the T = 
contribution, respectively, and the second terms yield the 
finite temperature corrections. Both yield the universal 
critical contribution to the free energy. T CI . Note that 
upon rescaling r — > 7'A 2 , T — > TA Z and ry Q — > r\ a A Zo, ~ z '', 
the dependence on the cutoff A reduces to a global multi- 
plication by A d+Z with d = 2, that can be later absorbed 
into the trivial scale dependence of the free energy. Fur- 
thermore, the 77 dependence of the universal kernels re- 
duces to 

/oo,«(T, r, Va ) = —fZ%(^T, r), (C9) 

which simply follows after substitution u> — > uj/rj a in the 
integral of Eqs. ([C7j and ([C8]) . 

We distinguish between a longitudinal and transversal 
part of the free energy, JF cr = J2 a -^cr,a, with 



= / dA e -^ A ^-/- d Q (r, Q ( e A )T( e A ),r(e A )). 

JO Va\ e ) 

(CIO) 

The running dynamical parameters r] a (b), temperature, 
T(6), and control parameter, r(b), entering Eq. (|C10[) are 



15 



given by the RG equations 



9?b 



d log b 



(z> - z)»7>> 



drj. 



d log b 
dT 



~— = ( Z< - z)t 1< 



= zT, 



(Cll) 
(C12) 
(C13) 



d log b 

f) r o 1 f) 

«£{<,>} 



du 
d log 6 



= (4 — d — ,z)u 



327rr/< 



(C15) 



where d = 2, z< = 2. and z> = 3. The initial condi- 
tions are rj a {l) = r] a , T(l) = T, r(l) = r, and u(l) = u. 
The temperature dependence of the RG flow of the quar- 
tic coupling, u, can be neglected in the asymptotic low- 
energy limit so that Eq. (|C15|) coincides with Eq. (|28bl) . 

One can convince oneself that, of course, the critical 
properties of physical observables do not depend on the 
specific form of frequency scaling employed in the RG 
process, i.e., on the value of the exponent z. In par- 
ticular, note that the combinations of scaling variables 
u(b)/rj a (b) and r) a (b)T(b), arc independent of z. As a 
consequence, the flow of the running control parameter 
r(b) is also independent of z as it is only determined by 
these combinations. Moreover, the additional multiplica- 
tive factor l/-q a (e x ) in the integrand of Eq. (|C10|) just 
cancels the dependence of the trivial scaling dimension 
of the free energy on z so that J- cr is also independent of 
the specific choice of z. 

The correlation length is obtained from r(b) with the 
help of the identification £~ 2 = r(b)/b 2 \ b _ A ^. The zero- 
temperature part of the transversal integration kernel, 
Eq. (|C7|) . proportional to r 2 governs the T = How of 
the control parameter. Only this leading term was kept 
in Eq. (j2"5a|) . The neglected longitudinal kernel, Eq. (|C8|) . 
gives a subleading nonanalytic correction to £ -2 on the 
order 0(u e er 3 ^ 2 ) at T = with the quartic coupling 
u e ft/ri < = m(A£)/t7<(A£) at the scale b = A£. On the 
other hand, the correlation length at finite temperatures 
in the regimes (II) and (III) of the phase diagram in Fig. [2] 
is controlled by the extended quantum-to-classical cross- 
over. This crossover is associated with scales b obey- 
ing the relation, 7?<(6/A)r(6/A) < 1 < ry>(6/A)T(6/A), 
which corresponds to an extended scaling range due to 



the presence of different dynamical exponents, z > =/= z<. 
Solving for £ in these regimes, one recovers the self- 
consistent Eq. (|3"T|) . 

The critical thermodynamics presented in Sec. IIII CI 
follows from expression (|C10p for the free energy. In par- 
ticular, the critical compressibility in regimes (I) and (II) , 
see Eq. (|4T)) . derives from the T = part of the transver- 
sal integration kernel [Eq. (|C7I) ]. 



1. Scaling variables rj > and n < 



The RG Eqs. (jCTT]) and JUTJ) imply that either 1 /?/> 
or r/ < is irrelevant when fixing the scaling dimension of 
the other dynamical parameter to zero, z = z < or z = z>, 
respectively. Here, we would like to demonstrate that 
these parameters are in fact dangerously irrelevant. For 
this purpose, it is sufficient to show^i that the critical 
thermodynamics is not well-defined in the two limits, 
l/?7> — > and ?;< — > 0. First, consider the contribution 
to the critical free energy deriving from the perpendicular 
mode. At T — 0, it simplifies to 



1 



cr,< lr = o 



32tt?7 



< Jo 



d\e~^[r(e A )} 2 . (C16) 



From this expression follows the aforementioned critical 
compressibility in Eq. (|T7|) . As JF CTi < is proportional to 
the inverse of r/ < , critical thermodynamics is obviously 
not well-defined in the limit r/ < — > 0. Second, consider 
•^cr > deriving from the longitudinal mode at criticality 
r = 0. The leading temperature dependence is given by 



V> Jo 



(C17) 



The dependence on r/> and T obtains after the appropri- 
ate substitution A — > A — ^ log(ry>T) in the integral and 
taking the limit of infinite cutoff, A — * oo, which elimi- 
nates the residual T dependence due to the lower bound- 
ary of the integration range. From Eq. (|C17[) follows the 
leading behavior of the specific heat in regimes (II) and 
(III), Eq. @2J). Clearly, it follows from Eq. ^Tt\ that 
thermodynamics is also ill-defined in the limit ?/> — > oo. 
Taken together, this identifies i] < and l/?7> as danger- 
ously irrelevant variables for the respective choices of z. 
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